Quantum transport efficiency and Fourier's law 
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We analyze the steady-state energy transfer in a chain of coupled two-level systems connecting two 
thermal reservoirs. Through an analytic treatment we find that the energy current is independent 
of the system size, hence violating Fourier's law of heat conduction. The classical diffusive behavior 
in Fourier's law of heat conduction can be recovered by introducing decoherence to the quantum 
systems constituting the chain. Implications of these results on energy transfer in biological light 
harvesting systems, and the role of quantum coherences and entanglement are discussed. 
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In recent years, energy propagation in systems that 
must be described in a quantum mechanical way has be- 
come a growing field. This growth is partially due to the 
fact that the understanding of how energy flow can be 
controlled and efficiently distributed has been identified 
as one of the crucial fields of study for the development 
of modern societies [2 [2] . One of the conceptual pillars 
in energy transport, the validity of Fourier's law of heat 
conduction, has become an active area of investigation 
and has been investigated in classical |3l |4j and quantum 
systems [5H7]. 

Since experimental evidence for quantum coherent ex- 
citation transport in the early light-harvesting step of 
photosynthesis has been presented [H [S] , investigations 
in systems of molecular biology have focused on the ques- 
tion, to what extent quantum mechanics contributes to 
the near perfect transport efficiency in light-harvesting. 
The emphasis has been put on the transient transport 
efficiency of an initial excitation in the presence of noise 
and disorder [I0l - [T2] . The experiments have been per- 
formed with pulsed femtosecond laser sources to excite 
and probe the molecule samples, whereas it has been 
suggested [I4HI6] that the light-harvesting process in 
vivo would be described more accurately in a steady- 
state scenario, because the light flux coming from the 
sun is essentially static on time scales that are relevant 
for molecular excitation transport. Here, we adopt the 
steady-state view on excitation transport in photosyn- 
thetic light-harvesting complexes, and evaluate the im- 
pact of Fourier's law on transport in quantum systems of 
this particular kind. 

An important step in the understanding how Fourier's 
law emerges from the quantum domain has been done 
by Michel et al. [6]. In this work Fourier's law is de- 
rived for a model system that is a chain of N identical 
coupled subunits, where each of the subunits has a sin- 
gle ground state and a narrow "band" of equally spaced 



FIG. 1. Chain of two-level quantum systems with its terminal 
sites coupled to heat baths of different temperature. 



excited states. In the present work, we employ a simi- 
lar system, i.e. a one-dimensional chain of two-level sys- 
tems, for which we compare the energy current in the 
classical analogue, where Fourier's law applies, with the 
quantum case, where we find the energy current to be 
independent of the chain length. This means that for the 
one-dimensional chain of two-level atoms Fourier's law 
applies for the classical variant but there is a distinct 
violation in the quantum transport scenario. By intro- 
ducing dephasing to the quantum model, we can study 
the transition from coherent to incoherent transport and 
show how Fourier's law can be recovered from the quan- 
tum case. 

Fourier's law of heat conduction states that the heat 
current through a classical macroscopic object is propor- 
tional to the applied temperature gradient |17j . 



J = -kVT, 



(1) 



where k, is the thermal conductivity. For a one- 
dimensional homogeneous object, the heat current is 
therefore determined by the temperature difference of the 
two heat baths AT, and the object length L. Generally, 
the validity of Fourier's law does not seem to be strictly 
linked to the classical or quantum nature of the system. 
For example, in the classical limit, for diffusive systems 
Fourier's law can be applied, but for ballistic systems in 
one and two dimensions there are divergences of the ther- 
mal conductivity as k ^ (see ^ for a review of heat 
transfer in low dimensional systems). For a discretized 
object composed of N equally spaced parts (sites), L (x N 
and thus 
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where c is a constant of proportionahty. For some one- 
dimensional quantum systems, one the other hand, there 
is evidence that Fourier's law is valid, i.e. a = [HIS]. 

The quantum system considered is a one-dimensional 
chain of > 2 two-level systems with coherent next- 
neighbor couplings as depicted in Fig. [T] The Hamilto- 
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where ct|, CTj!", and cr^ are the Pauli-z, raising, and low- 
ering operators in the basis of ground and excited state 
of the /cth two-level system, respectively, with on-site en- 
ergy huj and coupling strength g. Similar simple mod- 
els of coupled effective two-level systems are used in re- 
cent analyses of energy transfer in photosynthetic com- 
plexes [TDHT^ . The influence of the two heat baths is 
modeled by incoherently coupling each of the terminal 
sites to a bosonic heat bath described by a master equa- 
tion of Lindblad form. The system dynamics is then de- 
scribed by the master equation 



(4) 



where Ck acts on the first (last) site for k = 1{N), re- 
spectively, and is given by 

CkP = Vk{nk + l) {^kP^k - \ {(^k ^k^P]^ 



The first term in Ck accounts for emission into the reser- 
voir, the second term accounts for absorption, F^: is the 
interaction rate, and Uk = 1/ [e:iip{huj / {ksTk)) — 1] is the 
temperature-dependent mean excitation number at the 
resonance frequency in the respective bosonic thermal 
reservoir |13) . with fc^ being Boltzmann's constant. 

The expression of the heat current for a quantum sys- 
tem, Jq , is derived from the time-derivative of the energy 
of the system. 
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which vanishes in the steady state. When inserting Q 
into this expression, we obtain 



= Tr {HCip + HCnp) =■ Ji + Jw, 



(7) 



on the basis of which one can define the heat current 
to/from the respective reservoirs, both being of oppo- 
site sign, but equal in magnitude [T3|. The heat current 
through the chain is therefore equal to the net energy that 
enters the network from one reservoir and exits to the 
other per unit time, i.e. the quantity Jq = \Ji\ = \ Jn\- 



A straightforward evaluation of Jq for our system in the 
steady state yields the compact expression 

Jq = 7ifia;(si - (cTj^crf )) - ^^^(((Ti^cra") + (o-]"cr^)), 

(8) 

where 71 — ri(2rti -f 1) denotes the effective coupling 
to the reservoir, si = ni/{2ni + 1) is the excited-state 
population of a single two-level system in thermal equi- 
librium with reservoir 1, and all expectation values are 
taken with respect to the steady state of the chain. The 
heat current in the steady state is thus solely character- 
ized by the excited-state population of the first site and 
its specific energy gap, and since — {'^1^2)*^ 

is furthermore given by the real part of the coherence 
between sites one and two. An analogous expression can 
be given for the last site of the chain, which is connected 
to the second heat bath. 

For the complete expression of the heat current, 
we need the excited-state population of the first site, 
(a^ai), and the coherences between the first two sites, 
(crj^cr^). The excited-state populations of the individual 
sites in the steady state can be obtained from consid- 
ering specific matrix elements of the master equation of 
the kind fticT-^a^) = Tr(cr+cr^p) = 0. There are differ- 
ent cases: sites 1 and N, which are connected to their 
respective heat bath, and the remaining sites, which are 
in the middle of the chain. The relevant equations for 
the terminal sites k = 1 and k = N yield: 

71 (si - {crtcTi)) = ig{{a+(72) - (o-f cr^)). 



(5) For the inner sites, 1 < fc < iV, we obtain 



(9) 



that is, the imaginary part of all coherences between 
neighboring sites are equal. These equations motivate 
the following general form for the excited-state popula- 
tions of the terminal sites: 



> si - A/71, (a+(7jv> ^SN + A/-fN. (10) 



The transport along the chain thus causes a shift of the 
excited-state population of the terminal sites from the 
thermal equilibrium by A/7fc, where A — ig{^{a^ a^) — 
(6) {^i^ty)- The coherences, and thereby A, can be ob- 



tained by a similar argument. Summing up coherences 

of the steady state, ^ Y^k=i ^^t ^k+i) = ^-i provides the 
equation 



-ig ((cTj^CTi ) - {cft^a^)) = ^(cT]^ ) + ^( 



IN , + 

^\'^Af-l^JV/' 



the imaginary part of which, using ([9| and ([10]), yields 
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Next, for the complete solution of the heat current, we 
need the real part of the next-neighbor coherences, i.e. 
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FIG. 2. (Color online) Heat current as a function of the sys- 
tem size in a log-log plot (left) for the classical model for 
different values of the coupling V, and (right) for the quan- 
tum model for different values of the dephasing rate 7. Un- 
less otherwise indicated V — 1, g — 1, ksTi = 1, huu = 1, 
Ti = Tiv = 1, and Tn = 0. 

{'^1^2) + {'^i '^t) ■ For a chain of sites with uniform 
on-site energy as considered here, these coherences are 
purely imaginary, which can be shown from the structure 
of the master equation (see appendix). We thus arrive 
at the final expression of the heat current for a uniform 
quantum chain 

Jq — fi.ajA = — 2fiii;gIm(cr^(T^). (12) 

Here, we observe an important property of the heat 
current. It is independent of the chain length and thus 
violates Fourier's law, that is, the thermal conductivity 
scales as K ~ A^, and thus a = 1. 

If one considers technical applications, this result may 
open a new perspective to design materials with a more 
efficient heat transfer, free from losses due to the sys- 
tem size. Furthermore, provided that energy transport 
in light-harvesting complexes is governed by steady-state 
properties, this fact also provides a novel perspective on 
the understanding of highly efficient energy transport. 

Next, for comparison of the quantum model with the 
analogous classical model, we derive the heat current for 
the latter, which corresponds to the symmetric simple 
exclusion process jiTS], or Forster-type hopping [19\ It is 
a chain of N sites, which each may carry a single particle 
(excitation) that probabilistically moves between neigh- 
boring sites. This diffusive model fulfills Fourier's law. 
The classical probability for a particle to be at site k 
is given by Pk- The master equation ^ is thus turned 
into a Pauli master equation, i.e. a set of classical rate 
equations: 

Pi ^Fim + Pi{- Tiim + 1) - Tim -V)+ VP2 

Pk =V {Pk-i + Pfe+i - 2Pfe) {k^l,N) 

Pn ^^NriN + Pn{ - T^N{nN + 1) - ^NnN - V) + VPn-i 

where V is the constant rate to hop between sites. 
In the classical system, the heat current is defined by 
Jc — \y {Pi+i — Pi) \ 7 i-e. the net transfer rate of energy 
between sites, which in the steady state yields 

_ lijNVjsi ~ sn) 

V{ji + jn) + 1i1n{N ~ 1) ■ 
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FIG. 3. (Color online) Heat current as a function of tem- 
perature Ti with Tn = fixed, for the classical transport for 
different sizes of the system (left), and the quantum transport 
for = 4 and different dephasing ratios (right). Parameters 
are fcs — 1, fko = 1, Fi = Fjv = 1, V = 1, and g = 1. 

In the limit N oo, the heat current scales with the 
system size as Jc ~ V{si — SAr)/iV. Therefore, the heat 
current of the classical analogue obeys Fourier's law with 
K = const, and a = 0, in contrast to the quantum system. 

The difference between the heat current of quantum 
and classical systems can be lifted by adding a dephasing 
environment to each site of the quantum model. This 
amounts to introducing an additional term for every site 
in the master equation (|4]): 

k=l ^ ^ 

The results for the classical and quantum chain are given 
in fig. [2] on a log-log scale. The classical model features 
a linear dependence in the system size, for high enough 
values of N as expected. The heat current of the quan- 
tum case without dephasing is also linear in iV, but con- 
stant. However, when additional dephasing is applied, 
the heat current is suppressed and now features a size- 
dependence as for sufficiently large values of the 
dephasing rate 7, as confirmed by a numerical analysis 
of fitting the heat current to a power-law (see appendix) . 
By adding dephasing to the quantum system, we can thus 
recover the classical 1/A^-dependence of the heat current. 

In a common interpretation of a dephasing environ- 
ment, dephasing is caused by fluctuations of the on-site 
energy of every site. The excited state then effectively 
forms a band of states that is separated by a gap from the 
ground state. Adding dephasing thus effectively recovers 
the quantum model treated in [B], and yields the same 
qualitative result concerning the validity of Fourier's law 
regarding its dependence on the system size. 

Turning from the system size to the temperature- 
dependence, we find that in the quantum system the heat 
current features a strong dependence on the tempera- 
tures of the heat baths. Fig. |3] collects the temperature- 
dependencies for both models. The heat current of the 
classical system saturates for high values of the tempera- 
ture, which constitutes a violation of Fourier's law. This 
is due to the fact that the system has only two levels, 
which implies a finite heat capacity of the system. There- 
fore, it cannot transport an arbitrary large amount of 



4 




FIG. 4. Region of parameter space where the steady state 
may exhibit entanglement for 71 = 7jv. Within the shaded 
area of all possible values for Si and sjv, the upper (darker) 
parameter region indicates values for si and sjv where entan- 
glement can occur (dashed boundary not included). That is, 
only for si and S]v in this region, parameters g and 71 = 7jv 
exist such that the steady state is entangled. Entanglement 
cannot occur for any values of the coupling parameters out- 
side the darker shaded region. 



energy, and thus cannot scale linearly with the temper- 
ature for a large temperature difference. The quantum 
transport features a more intricate behavior. For a high 
temperature of the hot heat bath, its mean number of ex- 
citations and thereby 71 increases, causing a Zeno-type 
effect that reduces the transport efficiency of the system. 
With additional dephasing, the temperature-dependence 
of the heat current of the quantum system approaches a 
qualitatively similar saturating behavior as in the classi- 
cal system. 

Considering excitation transfer processes in light- 
harvesting complexes of photosynthesis, there are several 
lessons to be drawn from our model. Assume that, in- 
dividual sites of the chain represent the pigments that 
can carry an excitation, e.g. chlorophyl molecules, and 
the hot and cold reservoirs, which supply and consume 
excitation, correspond to the radiation field and pho- 
tosynthetic reaction centers, respectively. Even though 
the coupling topology of the network of pigments is not 
necessarily linear, we can arrive at a number of princi- 
ple observations: There exist quantum coherences in the 
system in the steady state that are static. It is these 
coherences alone that determine the transport properties 
as indicated in (12). Secondly, the violation of Fourier's 



law for the quantum transport provides a new way of 
how Nature may exploit quantum mechanics to gain size- 
independence, and thus for larger complexes to achieve a 
more efficient excitation throughput. 

A relevant point in this respect is the influence of dis- 
order in the system, and the observation that additional 
noise may unlock the effect of localization in disordered 
systems for transient transport processes [101 E] ■ How- 
ever, figs. [2] and [3] show that additional noise due to local 
dephasing reduces the observed heat current. Although, 
here, this result is obtained for a chain with uniform on- 
site energies and inter-site couplings, i.e. in the absence 



of disorder, we have also numerically investigated disor- 
dered chains. To this extent we have sampled the heat 
current in chains with N — 5, with all on-site energies 
hujk and couplings gki randomly chosen from a uniform 
distribution in the interval [0, 1]. In 8662 of the 10 000 
disorder samples, we found dephasing to reduce the heat 
current. Whenever additional dephasing is found to in- 
crease the heat current, the original random configura- 
tion exhibited a heat current below average of the entire 
random ensemble. We thereby extended what has been 
observed in the transient case [T^] to the steady-state 
scenario. 

With the perspective of identifying conceivable biolog- 
ical realizations of this transport scenario, an interest- 
ing aspect is the question whether entanglement is gen- 
erated and what role it plays, as addressed in [201 
Fig. [4] summarizes for which parameters entanglement of 
the non-equilibrium steady state occurs for a chain with 
N = 2 and equal effective bath rates 71 = jn- We find 
that entanglement can occur, but only in specific regions 
of the parameter space. Furthermore, for rates Fi = F^v, 
the steady state is never entangled for any choice of bath 
temperatures and coupling g. A bias in the bath rates, 
however, may drive the system to an entangled steady 
state. Depending on the interaction strength between 
the sites, entanglement may exist for certain range of 
temperatures. In contrast to entanglement studies in 
photosynthesis in the present scenario we find that 
the occurrence of entanglement is not equivalent to, and 
does not necessarily come with the mere presence of co- 
herences. It is thus an additional feature. 

To conclude, we have analyzed the energy transfer in 
a quantum system, formed by a paradigmatic chain of 
two-level systems, for which we found the heat current in 
the steady state to be independent of the chain length. 
This finding constitutes an important addition and gives 
a new twist to the study of the validity and emergence 
of Fourier's law in classical and quantum systems. We 
recover Fourier's law in the quantuni-to-classical transi- 
tion by adding dephasing that destroys quantum coher- 
ences. On the other hand, the violation of Fourier's law in 
the considered system encourages applications of quan- 
tum systems for efficient energy transport in specially 
designed materials by exploiting the size-independence. 
From the viewpoint of (excitation) energy transport in 
biological systems, the violation of Fourier's law consti- 
tutes another possibility of how quantum effects may be 
exploited in light-harvesting systems. It is the coherences 
in the system that govern the transport properties by de- 
sign, whereas entanglement may appear independently 
and in addition for a sufficiently large non-equilibrium. 
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I. APPENDIX 

A. Purely imaginary coherences 

From the structure of the master equation in Lind- 
blad form (4), an ordinary linear differential equation, 
one can directly infer that next-neighbor coherences are 
imaginary. In Liouville space the equation reads p = Lp, 
where p is the vector of all matrix elements, which are 
coupled linearly by the matrix L, the Liouvillian. For 
our purposes, it is helpful to introduce notation for the 
matrix elements: 

12. ..N{ik . . . q\p\jl ■ ■ ■ r)i2...N = Pij,ki,---,qr, 

where indices are grouped by subsystem. Vectors are 
products of basis vectors of the individual sites with 
ground state |0) and excited state We thus treat 
matrix elements with possible indices "0" and "1" . 

The LiouvUian L is a sum of three parts, each of which 
couples certain matrix elements, which yields indepen- 
dent sets of coupled matrix elements. It is possible to 
distinguish independent sets by observing general rules 
of how the Liouvillian couples matrix elements. We for- 
mulate these rules by the way indices are transformed by 
the Liouvillian. 

The Lindblad terms Ck of the Liouvillian inject or 
extract excitations at the terminal sites of the chain. 
Thereby, they transform matrix elements into one an- 
other that differ only by a pair of "00" and "11" indices 
of the first/last subsystem, e.g. poo,or ^ Pii.or- This 
constitutes a change of the total number of indices "0" 
and "1" by two, hence leaving the respective total number 
of indices "0" and "1" even or odd. Since H commutes 
with the excitation number operator, the commutator 
that appears in the Liouvillian leaves the total number 
of excitations invariant and hence couples only matrix el- 
ements with the same number of indices "0" and "1" , re- 
spectively. The coherent dynamics captures the exchange 
of excitations between neighboring sites and, in terms of 
matrix elements, couples those that can be transformed 
into each other by exchanging a "0" and a "l"-index be- 
tween neighbors, while maintaing the relative index po- 
sition, i.e. left and right indices are transformed within 
themselves, e.g. por.io ^ Pii.oo- This implies that the 
ground state poo.oo,... is coupled to all other populations, 
i.e. matrix elements with indices of the form pu.jj^..., and 
only to those coherences that contain an equal number of 
indices "1" on the left and right. The remaining matrix 
elements form an independent closed set of equations, 
whose steady-state solution is therefore the trivial solu- 
tion, where all matrix elements vanish. (The set of equa- 
tions that includes the populations is not solved by the 
trivial solution in the steady-state because it is subject 
to the boundary condition Trp = 1.) Note, that the di- 
agonal of L contains only coefficients with negative real 
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parts, meaning that all matrix elements would decay to 
zero if not sufficiently maintained by a positive contribu- 
tion due to another element. A population is coupled to a 
next-neighbor coherence, e.g. pn^o ^ Poi.iO) with a cou- 
pling ±ig such that (real) populations pump the imagi- 
nary part of the next-neighbor coherences, and vice versa. 
The contributions to/from the real part of the latter can- 
cel, leading to their decay. In longer chains {N > 2) next- 
neighbor coherences are also coupled to next-to-ncarest 
neighbor coherences, e.g. poi, 10,00 ^ Poi.oo.iOi with the 
same factor ±ig thus coupling imaginary (real) part of 
the former to real (imaginary) part of the latter. There- 
fore, imaginary and real parts of next-neighbor coher- 
ences belong again to different and independent sets of 
coupled differential equations. The real parts have the 
trivial solution, whereas the imaginary part is non-zero 



in the steady state. 

B. Fit of dephasing numerics 

To analyze the behavior of the quantum system under 
the effect of dephasing, we take the logarithm of (2) and 
obtain a linear dependence between log J and log N. 

log J = log (cAT) + {a-l) log N. 

By a linear regression of the numerical data of fig. 2, for 
7 = 5 we obtain a value a = 0.0242 with a regression 
coefficient R = 0.999985f. The small discrepancy with 
Fourier's law {a = 0) is due to the finite-size of the sys- 
tem. 



